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This Letter describes the search for a new heavy charged gauge boson W' decaying into an electron 
and a neutrino. The data were collected with the DO detector at the Fermilab Tevatron pp Collider at 
y^s = 1.96 TeV, and correspond to an integrated luminosity of about 1 fb _1 . Lacking any significant 
excess in the data in comparison with known processes, an upper limit is set on a w i x B(W' — > ev), 
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and a W' boson with mass below 1.00 TeV can be excluded at the 95% C.L., assuming standard- 
model-like couplings to fermions. This result significantly improves upon previous limits, and is the 
most stringent to date. 

PACS numbers: 12.60.Cn, 13.85.Rm, 14.70.Pw 



The standard model (SM) describes the fundamental 
fermions and their interactions via gauge bosons at a high 
level of accuracy, but it is not considered to be a complete 
theory. Additional gauge bosons are introduced in e. g. 
left-right-symmetric models (broken SU(2)l x SU(2)r) 
or in grand unified theories which may also involve su- 
persymmetry (e.g. Eq) Assuming the most general 
case, a new gauge group can comprise a new mixing an- 
gle £, new couplings to the fermions g' and a new CKM 
matrix U' . In some models the W boson (W' + or W'~) 
is right-handed, and decays therefore into a right-handed 
neutrino and a charged lepton. However, such a neutrino 
has not yet been observed. The mass limits on the W 
boson will in general depend on £, g', U' , and the masses 
of possible additional neutrinos. 

In this Letter we make the assumption that there is no 
mixing, g 1 is equal to the SM coupling, U' is equal to the 
SM CKM matrix, and that the decay channel W -> WZ 
is suppressed. Furthermore, the width Tw' of the W' 
boson is assumed to scale with its mass myyi , 
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The factor of 4/3 is applied in order to account for the 
decay into the third quark family (e. g. W — > tb) which 
is possible for mw' above the kinematic threshold for this 
process. In case of the existence of additional generations 
of fermions, it is assumed that they are too heavy to be 
produced by a W decay. This generic model has been 
introduced by Altarelli et al. 0]. It corresponds to the 
manifest left-right symmetric model [|| with light right- 
handed neutrinos if the W boson is right-handed. In this 
report, the general approach Q is considered, where the 
additional gauge boson W can be right- or left-handed. 

The W boson has been searched for previously by the 
DO HHEl] and the CDF experiments in various 

final states. The most restrictive limit so far is mw > 
800 GeV at the 95% C.L. reported by DO (W -> qq 1 , 
Run I). 

Data collected with the DO detector [l(| at the Fer- 
milab Tevatron pp Collider at a center-of-mass energy of 
1.96 TeV are analyzed for the production of W' bosons 
and the subsequent decay into an electron and a neu- 
trino. The neutrino can not be detected, but it gives rise 
to missing transverse energy ($t) m t ne detector. The 
dataset corresponds to an integrated luminosity [TJJ of 
0.99 ± 0.06 fb -1 , and was collected between 2002 and 
2006 during Run II of the Tevatron. 

The DO detector has a central-tracking system, con- 
sisting of a silicon microstrip tracker (SMT) and a cen- 



tral fiber tracker (CFT), both located within a 2 T su- 
perconducting solenoidal magnet, with designs optimized 
for tracking and vertexing at \g\ < 3 and \g\ < 2.5, re- 
spectively, where i] = — In tan | is the pseudorapidity 
and 9 the polar angle w.r.t. the z-axis (proton-beam di- 
rection). A liquid-argon and uranium calorimeter has a 
central section (CC) covering pseudorapiditics \r)\ up to 
w 1.1, and two end calorimeters (EC) that extend cov- 
erage to |?/| w 4.2, each housed in separate cryostats. 
An outer muon system, at \r/\ < 2, consists of a layer 
of tracking detectors and scintillation trigger counters in 
front of 1.8 T iron toroids, followed by two similar layers 
after the toroids. Luminosity is measured using plastic 
scintillator arrays placed in front of the EC cryostats. 

Different SM processes contribute to the electron and 
Et final state: inclusive production of W or Z bosons, di- 
bosons {WW, WZ, ZZ) or it pairs where at least one bo- 
son or one top quark decays into electrons directly or via 
tau decays. In these processes the missing energy is due 
to the neutrino. There are also two sources of misiden- 
tification background that can contribute to the electron 
and |?t final state: QCD multijet background with one 
jet misidcntificd as an electron and energy mismcasurc- 
ment which can cause large $t cither along or in the op- 
posite direction of the electron, and Z — > ee events where 
one electron is lost (e. g. entering non-instrumented sec- 
tions of the calorimeter) or misreconstructed. The latter 
case can lead to large Et- 

The W' signal and SM processes (including Z — > ee) 
have been simulated with the pythia 6.323 [l2[ Monte 
Carlo program using the CTEQ6L1 [lj| parton distri- 
bution functions (PDFs), except for the QCD multijet 
background, which is estimated from data. The gener- 
ated events are passed through a detailed detector simu- 
lation based on GEANT 14| , and combined with randomly 
triggered events from data to simulate the effects of pile- 
up and multiple interactions. Higher order corrections to 
the pythia leading order cross sections (K factors) have 
been applied. The next-to- next-to-leading order (NNLO) 
K factors and errors due to PDF uncertainties for the sig- 
nal, the W and the Z samples are extracted from Ref. 
0; the NNLO (NLO) cross section for tt (di-boson) 
production is taken from Ref. [l6[ ([I3])- 

The signal cross section falls steeply with increasing 
mass of the W' boson. In addition, for very large masses 
the on-mass-shell production of W' bosons is heavily sup- 
pressed due to the smallness of the PDFs at large x. As 
shown in Fig. [l] the Jacobian distribution no longer ex- 
hibits a pronounced peak. The transverse mass mj is 
calculated from the transverse energy of the electron, £?fJ, 
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FIG. 1: Transverse mass mr distributions for different masses 
of the W boson (generator level, pythia). The event num- 
bers correspond to an integrated luminosity of 1 fb _1 . 



the missing transverse energy, an d the azimuth angle 
[3] difference between the electron and Et vi a 



rriT 



2E^E T {1 - cos A</>(electron, E T )). (2) 



Events triggered by a set of inclusive single electron 
triggers are considered. Electrons with > 30 GeV 
passing the offline identification criteria are selected. 
Monte Carlo studies have shown that the majority (« 
80%) of the electrons stemming from the W decays 
are emitted into the central detector region. Since 
the forward detector region exhibits a small signal-to- 
background ratio, only electrons reconstructed in the CC 
are used in the analysis. Electromagnetic clusters are 
built around a calorimeter seed. Such clusters consist of 
cells in a cone (AR — -J (Ary) 2 + (A0) 2 < 0.4) around 
the seed. Furthermore, the electron shower is required 
to be isolated in the calorimeter, and to deposit most 
of its energy (> 90%) in the electromagnetic part of the 
calorimeter. The isolation I = [E$£ — Bem]/-^emi which 
uses the total shower energy, E$£, in a cone of radius 
R = 0.4 and the electromagnetic energy, E^, in a cone 
of radius R = 0.2, is required to be less than 0.2. A cut on 
the electron shower shape variable is applied to separate 
electromagnetic from hadronic showers. The electron is 
required to have a track matched in z and 4> direction 
and to stem from the primary vertex. Correction fac- 
tors are applied in order to take differences in the recon- 
struction efficiencies observed in data and Monte Carlo 
into account. Finally, the energy dependence of the ba- 
sic electron reconstruction criteria has been studied with 
simulated electrons from W decays. The reconstruction 
efficiency is found to be constant (94 ± 1 %) and does 
not exhibit a visible energy dependence within the statis- 
tical uncertainties of the Monte Carlo samples. The $t 
is calculated from all calorimeter cells. Corrections are 
applied to account for the electromagnetic and jet energy 



scales. We require $t > 30 GeV. 

Since the transverse momentum of the neutrino is ex- 
pected to be balanced by the electron transverse energy 
in signal events, a selection on the ratio of the energies is 
applied, 0.6 < E^/Et < 1.4. This requirement reduces 
instrumental backgrounds from misidentified $t- Jets 
are reconstructed with the iterative mid-point cone algo- 
rithm (R = 0.5) [HI. If any jets with p T > 15 GeV are 
present in the event, we require A</>(jet, electron) < 2.8 
and A0(jet, $t) < 2.8. These selections remove events 
from QCD multijet production. 

The contribution from QCD multijet events is esti- 
mated using a control sample derived from data with 
the same kinematic cuts. In this sample, the electron 
candidate fails the shower shape requirement. The re- 
sulting events are normalized to the data sample. The 
scale factor for the entire QCD multijet sample is ad- 
justed in the low reconstructed transverse mass region 
(tot < 30 GeV), which is dominated by QCD multi- 
jet background events, such that the sum of the PYTHIA 
Monte Carlo prediction and the QCD multijet sample de- 
scribes the data as shown in Fig. [2ja). The data are nor- 
malized to W boson production and decay in the ev mode 
using the W peak region (60 GeV < Tot < 140 GeV, 
as shown in Fig. [2ja)) because many efficiency and ac- 
ceptance errors largely cancel in this ratio. We use the 
theoretical prediction for the W boson production cross 
section a w x B(W -> ev) = 2583±^ pb from Ref. [H|. 

Jets may be present in conjunction with a W boson due 
to higher order QCD contributions. Since pythia does 
not properly describe the transverse momentum distribu- 
tion of the W boson in such processes, this spectrum is 
separately reweighted in events with one, two and three 
jets in order to match the distributions observed in the 
data. This correction affects 10% of the W Monte Carlo 
events. The sample defined by the selection cuts men- 
tioned above contains 452,984 data events compared to 
454,000 ± 35,000 events expected from SM processes and 
instrumental backgrounds after applying all corrections. 

Two kinds of systematic uncertainties contribute in 
this analysis. The uncertainties of the normalization in 
the W peak region (4%), the cross sections of the SM 
processes (4-10%), the electron reconstruction efficiency 
corrections (2%) and the scale factor for the QCD mul- 
tijet sample (7%) affect only the global normalization. 
Uncertainties on the PDFs, electron energy scale and 
resolution, jet energy scale, decay width Tw of the W 
boson, and the reweighting of the transverse momentum 
of the W boson lead to changes of the shape of the dis- 
tributions. 

In order to study the effect of the electron energy scale 
and resolution, the electron energies have been varied 
within the known uncertainties. The variations of scale 
and resolution are performed independently. The $t 
is recalculated after varying the electron energy. The 
overall uncertainty on the event numbers is large for the 



6 




FIG. 2: Comparison between data and background prediction: (a) Distribution of the transverse mass mj; (b) distribution of 
the electron transverse energy Ex in events with mr > 140 GeV; (c) distribution of the ratio of electron transverse energy and 
Et in events with mr > 140 GeV. The signal is shown for two different masses of the W' boson. 



W sample (4%), but small for the W signal (< 1% for 
500 GeV < m w , < 1200 GeV). The uncertainty of the 
energy resolution is an order of magnitude smaller than 
the energy scale uncertainty. 

In order to study the PDF uncertainty, the Monte 
Carlo events which have been produced using CTEQ6L1 
PDFs are reweighted to CTEQ6.1M.xx (xx = 0, . . . , 40), 
making use of the CTEQ6.1M PDFs and the 40 error 
functions [3] . The overall uncertainty varies between 3% 
{mw = 500 GeV) and 8% {mw = 1200 GeV). For the 
W sample an uncertainty of 3% is derived. The width 
of the W boson is known to about 2% [2Qj . This can 
cause a shift (~ 4%) of the tail of the transverse mass 
distribution of the W boson. Finally, the jet energy scale 
has been varied, and the $t recalculated. The resulting 
uncertainty is below 1%. 

The tail of the spectrum (mj- > 140 GeV) is now con- 
sidered to search for W' — > ev. A good agreement be- 
tween data and background prediction can be observed 
as shown in Fig. [2fb, c). In Table [H the breakdown of 
the individual contributions of the various background 
processes is given, including expected numbers of signal 
events. Since we do not observe any significant excess in 
the data, an upper limit is set on the production cross 
section times branching fraction aw x B{W' — > ev). 

The limit is derived using a binned likelihood for 
the whole transverse mass spectrum with 140 GeV < 
mr < 1000 GeV. The individual shape-changing system- 
atic uncertainties (up and down variation) enter the limit 
calculation via individual histograms; bin correlations are 
taken into account. 

A Bayesian approach [2l[ is used to calculate upper 
limits on the cross section for different resonance masses. 
A Poisson distribution is assumed for the number of ex- 
pected events in each bin of the transverse mass distri- 
bution, as well as flat prior probabilities for the signal 
cross sections. The prior for the combined signal accep- 



TABLE I: Event numbers in the data compared to the back- 
ground prediction after applying the cut on the transverse 
mass m,T > 140 GeV. For the signal and background pro- 
cesses statistical and systematic uncertainties are given. 



Process 


Events 


stat 


sys 


Data 


967 






Sum of backgrounds 


959 


21 


90 


W -> ev 


875 


20 


90 


QCD multijet (from data) 


27 


2 


2 


other 


57 


3 


4 


W -* ev 








m w > = 500 GeV 


1169 


24 


86 


m w > = 600 GeV 


393 


8 


32 


m w , = 700 GeV 


147 


3 


13 


m w > = 800 GeV 


51 


1.1 


5.4 


m w , = 900 GeV 


19 


0.4 


2.4 


m w > = 1000 GeV 


7.4 


0.2 


1.1 


m w , = 1100 GeV 


3.4 


0.1 


0.5 


m w , = 1200 GeV 


1.7 


0.1 


0.2 



tance and background yields is a multivariate Gaussian 
with uncertainties and correlations described by the cor- 
responding covariance matrix. 

The observed and expected 95% C.L. limits on the 
production cross section times branching fraction aw' x 
B{W — > ev) are shown in Fig. [3] The lower bound of 
the theoretical cross section is used to obtain the mass 
limit. Hence, an additional heavy charged gauge boson 
with mass below 1.00 TeV is excluded at the 95% C.L. 

In summary, a search for a new heavy charged gauge 
boson W' decaying to an electron and a neutrino has 
been performed using 1 fb _1 of data collected with the 
DO detector in Run II. We do not observe an excess in 
the data, and set upper limits on the cross section times 
branching fraction, which are of the order of 10 - 40 fb 
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FIG. 3: The observed and expected 95% CL. limits on the 
cross section as a function of the mass of the W boson, in- 
cluding statistical and systematic uncertainties. The expected 
limit assumes a background-only hypothesis. The theoretical 
expectation is displayed with its uncertainty. Also shown is 
the DO Run I limit @|. 



for W boson masses of 500 GcV < m w > < 1200 GcV. 
Further, a lower limit on the mass of the W boson is de- 
rived, assuming that the new gauge boson as introduced 
in 0] has the same couplings to fermions as the SM W 
boson. We exclude a W boson with mw' < 1-00 TeV at 
the 95% CL. This result represents the most stringent 
limit on the mass of a charged heavy gauge boson beyond 
the standard model to date. 
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